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1 Introduction 

The mathematical description of compressible flows is important from the point of view of ap- 
plications, domains such as aerodynamics and geophysics are the most natural to be mentioned 
here. On the other hand, complexity of the equations describing the flow delivers very interesting 
mathematical challenges. In spite of active research in the field, we are still far from the com- 
plete mathematical understanding of compressible flows. The only general existence results are 
available for weak solutions with homogeneous boundary conditions [7], [12]. As far as regular 
solutions are concerned, we have so far only partial results assuming either some smallness of 
the data, or its special structure. The problems have been investigated mainly with homogeneous 
boundary conditions ([4], [17]). For the overview of the state of art in the theory one can consult 
the monograph [18]. 

From the point of view of the aforementioned applications it seems very important to inves- 
tigate the problems with large velocity vectors, which lead in a natural way to inhomogeneous 
boundary conditions. Due to the hyperbolic character of the continuity equation the density must 
be then prescribed on the inflow part of the boundary. Existence issues for such inflow problems 
are investigated in [9], [10], [19], [20], [21], [22] and [26]. The mentioned group of problems 
can be regarded as questions of stability of particular constant flows. 

In the present article we would like to examine the issue of stability of Poiseuille type flow in 
pipe-like domain in compressible regime. The unperturbed flow is a solution to the compressible 
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Navier-Stokes system with homogeneous slip boundary conditions for given constant friction, 
constant density and constant external force (gravitation-like term). The Poiseuille flow is a 
special symmetric solution to the incompressible Navier-Stokes equations in cylindrical domains. 
Here it is viewed as a solution to the compressible Navier-Stokes system with constant density 
and constant external force, parallel to axis of the cylinder, given by the pressure. Hence the 
pressure, unknown in the incompressible model, is recognized as a given force. Such change of 
'observer' looks acceptable from the mechanical point of view. Thanks to that interpretation we 
obtain a natural physically reasonable flow in compressible regime. The mathematical objective 
of this article is to establish stability of such flow under some structural assumptions limiting the 
magnitude of admissible perturbations. 

Let us define the system. We consider steady flow of a viscous, barotropic fluid in a bounded, 
cylindrical domain in IR 3 , described by the Navier - Stokes system supplied with inhomogeneous 
Navier slip boundary conditions. The complete system reads 

pv ■ Vi> — pAv — (p + v) Vdiv v + Vir(p) = p F in f2, 

div (pv) =0 in f2, 

n ■ T(v, 7r) • r k + fv ■ r k = b k , A; = 1,2 on T, (1.1) 

n ■ v = d on T, 

P = Pin on F in , 

where f2 = [0, L] x f2 with f2 C IR 2 of class C 2 , T denotes the boundary of f2 (see Fig.l), v is 
the velocity field of the fluid, p its density, p and v are viscosity constants satisfying p > and 
(v + 2p) > 0, / > is the friction coefficient which may be different on different components 
of the boundary T, it = n(p) is the pressure given as a function, at least C 1 , of the density and 
F is an external force. T denotes the Cauchy stress tensor of the form T(v,n) = 2pDv + 
z/divt>Id — 7iTd where D = | ( Vi> + Vi> T ) is the symmetric gradient. Next, n and r k are outer 
normal and tangent vectors to dil. Boundary data p in , b, d will be discussed later. The boundary 
T is naturally split into three parts: 

r = {x e dQ : v(x) ■ n(x) = 0}, 

T in = {x e dtt : v(x) ■ n(x) < 0}, (1.2) 
r o «t = {x G dQ : v(x) ■ n(x) > 0}. 

Thanks to the chosen geometry of the domain, the above decomposition is easily illustrated by 
the figure 1 . 

We shall say few words about the physical interpretation of the system (1.1), in particular 
about the choice of boundary conditions (1.1)3,4. We would like to model a flow through a pipe. 
We assume that the fluid obeys Navier slip conditions on the walls of the pipe (r component of 
the boundary), hence natural conditions on T are d = and b k = 0. However, the mathematical 
requirements impose a need to prescribe the boundary conditions on T in and Y out . From the 
physical viewpoint these parts are artificial, this is the area where the parameters of the velocity 
and density are measured. This gives us a freedom of choice of the type of boundary condi- 
tions on the inflow and outflow part, which can be fit to the mathematical approach, hence we 
choose inhomogeneous slip condition. Note that as the friction coefficient goes to infinity, then 
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Figure 1 : The domain 



the relations (1.1) 3 4 , at least formally, become the standard Dirichlet conditions describing the 
whole velocity vector at the boundary. Since the velocity does not vanish on the boundary, the 
hyperbolicity of the continuity equation impose a need to prescribe the density on the inflow 
part, which lead to the condition (1.1)5. The velocity field determines the characteristics of the 
continuity equation and in particular the total mass J n p dx is determined implicitly by (1.1)5. 
Our goal here is to analyze a perturbation of the Poiseuille type flow 



where x\ points the direction of the axis of the cylinder. It is one of the classical examples of lam- 
inar flows satisfying the incompressible Navier-Stokes equations in cylindrical domains. In the 
classical literature the flow is considered with homogeneous Dirichlet condition on the boundary. 
V p is then found as a solution to corresponding elliptic problem with Dirichlet boundary con- 
dition on each x\ - cut of Vt. Some explicit formulas on V p in certain domains are well-known 



In the case of slip boundary conditions that are subject of our analysis in this paper the 
flow (1.3) can be also found on each cut of the cylinder as a solution to elliptic problem with 
corresponding boundary condition (see Lemma 1 below). In certain domains it can also be 
expressed with explicit formulas (see [14] and the example below). Since we are interested in a 
general cylindrical domain, we will not have such formula but we show that the solution of the 
form (1.3) exists provided that Vt Q is sufficiently regular and, under the slip boundary conditions 
(1.1)3,4, V p does not vanish on the boundary. 

Lemma 1. Let VL^ = lx Q , where Qq C M 2 is smooth enough, fi > and f > 0. Then there 
exists a solution (V, Tl), such that V = (y p (x 2 , x 3 ), 0, 0), to the incompressible Navier-Stokes 
system with slip boundary conditions: 



v = \y p {x 2l x 3 ) 1 o,o], 



(1.3) 



([8],[11]). 




(1.4) 
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Moreover, there exist 6 = 9(f,fj) and a continuous function 00(f) such that 



V w = V P >6>0 in fioo. (1.5) 

and 

IIWIl. < (1.6) 
/i 

In addition if(V, Tl) is the Poiseuille solution, then (XV ', An), too. 

As we already said, the pressure in the Poiseuille flow can be regarded as an external force 
parallel to the axis of the cylinder. Like in the classical Poiseuille flow, we assume the pressure 
to be a linear function of x\. Hence it is natural to assume the form n = oo(f)x\, where 00(f) is 
a negative constant (the sign describes the directon of the flow). Then we see that V p , the first 
coordinate of V — (V p , 0, 0), is a solution to the elliptic problem 

/2AV P = U X1 = 00(f) < infi , 

H § £ + fV p = ondQo- ( } 

Now to prove Lemma 1 it is enough to apply the maximum principle to the system (1.7). We 
show the proof in the Appendix, at this stage we should have a closer look at the dependence 
00(f). This dependence can be justified considering the compatibility condition for the system 
(1.7), that reads J" 00 dx = f J dQo V p da. As we will see, the dependence 00(f) determines the 
assumptions we will have to make on the viscosity. 

Note that for 00 = the only solution is V p = 0, provided / 7^ 0. For / = (perfect slip) 
we should obtain a constant flow, hence we put 00 = for / = 0. The linear structure of (1.7) 
makes 00(f) continuous. Thus, we conclude 00(f) — > for / — > 0, and hence by (1.6) 

l|VV|| Loo ^0 for f^0. (1.8) 

This is an important conclusion since in order to show the energy estimate we have to control 
VV with the viscosity, and so we allow the viscosity to be low provided that the friction on T is 
small, what is a realistic assumption (see also the remarks after the formulation of Theorem 1). 
Let us illustrate the dependence 00 = 00(f) with the following example. 

Example. Take Q = B(0, 1) C IR 2 and ji = 1. Then, due to axial symmetry of the domain, 
it is natural to look for V p = v^(r) where r = \/x\ + x\. The boundary condition (1.7) 2 then 
reads v / + fv f \ r= i = 0. We require that Av = v rr + \v r depend only on /. Moreover, we expect 
to obtain a constant flow for / = (perfect slip) and classical Poiseuille profile for / = 00. The 
above considerations lead to the family of solutions 

where kf — ( 7r ~ 2 ^ +47r an d TF is the flux of the flow v- f through f2 . 
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For a perfect slip case / = we obtain a constant flow v = 2 ^y- and for a no-slip case 
/ = oo we get a classical Poiseuille profile v°° = TF(1 — r 2 ). On the boundary we have 
9(f) = vf(l) = what is a strictly positive constant. Finally, 



4/ f f 

u(f) = Av' = — — and Vir 



f + k, f + k, 

In particular u(f) < for / > and w(0) = 0. 

Before we formulate our main result, we need one observation concerning the boundary 
conditions. Note that, since V p is found on every x\ - cut of f2, we can impose the boundary 
conditions (1.4) 3; 4 only on r . On the other hand, in order to define small perturbations as a 
solution to (1.1) we have to measure the distance (in appropriate norms) between the solution to 
(1.1) and the Poiseuille flow V. Hence we need to consider the traces of the quantities from the 
boundary conditions of (1.1) with the function V instead of v. Since our analysis acts on a finite 
cylinder we define these traces at the bottoms T in U T out : 

6fc|r iB ur OTt = tir m ur out [n-T(V,f[)-T k + fV-T k ], d\ Fm = -trr in V p , d\ Tout = tr Tout V p , 

where tr denotes the trace operator. By (1.4), 

n ■ T(V, fl) • r k + fV ■ r k = and n ■ V = at T . 

Hence the construction of Poiseuille flow determines (k = 1,2): 

b k = d = on T . (1.9) 

The Poiseuille flow (V, IT) has constant density p, we set p = 1. Then in the chosen setting 
(V", IT) fulfills the following system 



pV ■ W - pAV -(p + v) Vdiv V + Vn(p) = -pw(/)ei in Q, 

div (pV) = in O, 

n ■ T(V,U) ■ r k + fV -r k = b k , A; = 1,2 on T, 

n -V = d on T. 



(1.10) 



We keep in mind that VII = oo(f)e\ and (1.9). 

We are now in a position to formulate our main result. To this end it is convenient to define 
the quantity which measures the distance of the data from the Poiseuille flow: 

A) = \\F + w(/)ei|| Lp (n) + d\\ w 2-i/ P ( T) + \\b k - b k \\ w i-i/ P ^ + \\p in - l||vin(r in )- (1-H) 

The main result of the paper reads 

Theorem 1. Assume that the boundary data of (1.1) is close to the Poiseuille flow (V, fl), more 
precisely, let Dq defined above be small enough. Assume that the friction f is large enough on 
T in andp > 3. Assume further that there exists a constant k > such that 

p > K,min{f\ ro , 1}. (1.12) 
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Then there exists a solution (v, p) G W£(Q) x Wp(Q) to the system (1.1) such that 

\\ v ~ V\\w*(n) + IIP - P\\w^n) < C(D ). (1.13) 
This solution is unique in the class of small perturbations of (V, p). 

Let us make some remarks concerning our main result. The condition on the viscosity (1.12) 
seems to be a serious constraint, but as we will see from the proofs we just need it to control the 
gradient of the Poiseuille flow, what yields this assumption natural (see also the remark in the 
proof of Lemma 6). We recall that VV" depends on the friction / on T , and in particular (1.8) 
holds. It follows that for small values of friction on r it is enough to assume that the viscosity 
is large enough, but only compared to the friction. This assumption is reflected in the condition 
(1.12). Theorem 1 admits the case of perfect slip /| r „ = 0, and in such case (1.12) reduces to 
p > 0, so no lower bound on the viscosity is required. In this case there is no bound on the 
size of V. However in this case V would be a constant flow. We shall recall that the friction at 
T in is chosen independently to / at r . From the point of view of modelling, the data at T in is 
given, hence it is important to focus the attention at r . The assumption p > 3 is required for 
the imbedding C [1], it is required to control the pointwise boundedness of Vt> and the 
density. 

Let us explain the main idea of the proof. We will follow an idea of Lagrangian type coor- 
dinates [5], [13], [16], [24]. The continuity equation is of hyperbolic type and contains a term 
u-Vw (where u and w are perturbations to the velocity and density introduced in the next section), 
which makes serious troubles for the issues of existence in case of inhomogeneous boundary con- 
ditions, see [9], [21], [18], [19], [20]. Here we overcome this obstacle by changing the system 
of coordinates in such a way that this term disappears (2.7). We obtain a more complex system 
but with structure suitable for an application of the Banach fixed point theorem. On the other 
hand our solutions are regular enough, thus we are able to go back to the original system keeping 
the well posedness of the original model. Our approach works since we are equipped with the 
maximal regularity estimate for a linearization of the equations in the Lagrangian coordinates - 
Theorem 2. This tool gives a complete control of the regularity of solutions. 

The rest of the paper is organized as follows. In Section 2 we introduce the perturbations 
as unknown variables obtaining the system (2.5). Next we introduce the Lagrangian-type coor- 
dinates that lead to the system (2.11) and we derive the necessary estimates for the Lagrangian 
transformation. In Section 3 we deal with the linearization of (2.11). For the linear system we 
show the estimate in Wp(fl) x Wp(Q). It is given by Theorem 2. The first step is the energy 
estimate (3.3). Then we consider the vorticity of the velocity and the Helmholtz decomposition 
to reduce the continuity equation to a sort of transport equation (3.22) that enables us to find the 
bound on || w\\ w^(n)- This result together with the properties of the Lame system lets us conclude 
Theorem 2. In the second part of this section we apply the estimates to solve the linear system 
and hence show that T given by (2.27) is well defined. In Section 4 we show the contraction 
principle for T. To this end we consider the system for the difference of two solutions and write 
it in a form (4.1) which has a structure of (3.1). The contraction results from the estimate (3.30) 
and bounds on the norms on of the r.h.s. of the system for the difference. At the end of Section 
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4 we apply the Banach fixed point theorem to solve the system (2.1 1) and conclude the proof of 
Theorem 1. 

Let us finish this introductory part with some remarks concerning notation. By C we shall 
denote a constant that is controlled, but not necessarily small. E shall denote a constant that 
can be arbitrarily small provided the data is small enough. Sometimes we will write E(-) to 
underline that we need the smallness of certain quantity. The functional spaces on f2 will be 
denoted without the symbol of the set, for example we will write Wjf instead of W£(Q) for 
standard Sobolev spaces of functions intergable with the p-th power with derivatives up to order 
k, wl~ l ^ p (dVt) denotes the Slobodeckij spaces, defining regularity of traces from Wp(Q), [1]. 
Finally, we will need to consider the density in the space Loo(0, L; L 2 (^o))- For simplicity we 
denote it as L 00 (L 2 ). We do not use different notation for scalar and vector valued functions, 
while matrix valued functions are written in bolded font. The coordinates of a vector are denoted 
by (•) > i.e. u = ( u W, u (2), u (3)). 

2 Preliminaries 

In this section we introduce perturbations of the Poiseuille flow (V,p) as unknown variables, 
what leads to the system (2.5). Then we introduce a change of variables that straightens the 
characteristics of the continuity equation. We obtain the system (2.11). The simplified form of 
the continuity equation in this Lagrangian framework makes it possible to apply the Banach fixed 
point theorem to the system (2.1 1). 

2.1 Reformulation of the problem 

We come back to the main system (1.1). Since we are interested in solutions that are small 
perturbations of (V", p), it is convenient to consider the perturbations as unknown functions. For 
technical reasons it is better to have u ■ n = on the boundary. Hence we start introducing 
Mo £ W£ (Q) such that u ■ n\r = d — d (recall that d = V - n).lt can be found as u = V0 where 
4> solves a Neumann problem. We assume that \\d — V ■ n\\ w i-i/ P ^ is small enough for 

v p + u^ ) \ n >e 1 (2.i) 

to hold for some 6\ > 0. In fact this is not really a restriction as we consider small perturbations 
of (V", p) and in Lemma 1 we have shown that V p is separated from zero. Now we take 

u — v — V — uq. (2.2) 

In particular we want our perturbed flow v to have the first component also separated from zero. 
This is quite natural constraint if we consider small perturbations of V. With the above definition 
of u this constraint reads 

V p + + 4 1} > 2 > 0. (2.3) 
Next we introduce the perturbation of the density (recall that p = 1): 

w = p-1. (2.4) 
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Substituting (2.2) and (2.4) to (1.1) and (1.10) we arrive at 

u ■ VV + V p d Xl u — fiAu — (fi + z/)Vdivw + 7 V w = F(u, w), 
V p d x , w + div u + (u + uo) ■ Vw = G(u,w), 

(2.5) 

n ■ 2/iD(w) • r k + f(u ■ T k )\ T = B k , 
n-u\ r = 0, w\r zn = w in , 

where 7 = 7r'(l) and 

F(u, w) = (u + M ) • V(u + u ) -uo-VV -V p d Xl u - [ir'(w + 1) - 7r'(l)] Vw 

-wu(f)e! + (w + l)(F + w(/)ei) - iu(u + m + V") • V(« + m + V"), 

w) = —w div u — (w + 1) div uo, 

B k = b k -n - 2fxD(V + u ) ■ r k - f(V + u Q ) ■ r k . 

From now on we focus on the system (2.5). Notice that V ■ VV = 0, hence the term — w(u + 
Uq + V) ■ Vm + u + V) is a higher order term and so the form of F and G implies immediately 
the following lemma. 

Lemma 2. Let F(u, w) and G(u, w) be given as above, then 

\\F(u,w)\\ Lp + \\G(u,w)\\ w i < C[(\\u\\ W 2 + \\w\\ w i) 2 + E(\\u\\ W 2 + \\w\\ w i) + D ], (2.6) 
where E denotes a small, compared to ji, positive constant. 



2.2 Change of variables 

With our smallness assumptions it is quite natural to solve (2.5) with a fixed point argument. 
However, a direct application of this method fails because of the nonlinear term u ■ Vw in the 
hyperbolic continuity equation. The idea to overcome this problem is to introduce a change of 
variables such that this awkward term vanishes. We look for the appropriate transformation as 
x = ip u+Uo (z) satisfying the identity 

V p d ei =V p d Xl + (u + Uo)-V x . (2.7) 

In the following lemma we construct the mapping ^« for arbitrary function u small in W% with 
vanishing normal component on the boundary r . 

Lemma 3. Let \\u\\ W 2 be small enough and u • n\ To = 0. Then there exists a diffeomorphism 
x = ipu(z) defined on Vt such that fl = i/)a(Q) and (2.7) holds with u + uq = u. 

Proof. A key point in the proof is the fact that V p > c > 0. In particular we are able to 
divide (2.7) by V p obtaining 
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where u — yp. Since u, V p G W 2 we have u G W£ and, since we are interested in small 
perturbations we can assume that 

\\u\\ W 2 « 1. (2.8) 

Now we can follow the proof from [20] and look for tp(z±,z 2 , 23) = ipz 2 ,z 3 (zi), where for each 
(z 2 , Z3) G r the function ip Z2 , Z3 is a solution to 

f 9«^£!*3 = 1, 9 a ^2!z 3 = « 2 (^*2,a3)> ^S!*3 = « (3) (^«2,«3), ( 2 9) 

I ^2,^( ) = (0,22,23)- 

Due to (2.8) we solve (2.9) for (z 2 , z s ) G Y in following [20] and show that there exists a set Q u 
such that ip(Qu) — > fl is a diffeomorphism. It remains to show that Q s = f2. To this end we 
examine the derivatives of ip. We have _D-?/> = Id + E where 



E = 





u 2 (ip(z)) d Z2 j* 1 u 2 (tp(s, z))ds d Z3 J 21 m 2 (^(s, z))ds 
u^^(z)) d Z2 tfuW^(s,z))ds d Z3 J^u^^(s,z))ds 



(2.10) 



and z := (z 2 , z 3 ). The first row of E reduces to since E^ = d Zi f* 1 ds = for i = 2, 3. Hence 

£^([1,0,0]) = [l,fi 2 (^)),«W(^))]. 

Now take rr n — > x G r and z n = 4>(x n ). Then we have 

lim n ^ OOJ DV(2n)([l,0,0]) = [l,u 2 (x ),« (3) (a*)]- 

The latter is parallel to T since tt • n|r = 0, hence we have <p(T ) = T (precisely we say it 
in a sense of tangent spaces since is defined only on f2). To examine the behavior of tangent 
vectors on r oMt notice that 

D^(z)([0,t 17 t 2 }) = ([0,f 1 (1>(z))M1>(*))]), 

where are given by appropriate entries of Dip, the important fact is that the first coordinate 
vanishes. Hence if we take x n — > x and z n = 4>(x n ), this time with x G V out , then 

lim n ^ oo Dip(z n )(l0,T 1 ,T 2 ]) = [0,f 1 (x ),f 2 (x )], 

what is parallel to T out . It shows that <p{Y out ) = T out . Since 0(T O ) = T by the definition of 0, 
we conclude that £l u = ^ and complete the proof. □ 

Now we proceed with transformation of the system (2.5). As a vector field u satisfying the 
assumptions of the Lemma we take u + uq where u is the solution of (2.5). So far we don't know 
if this solution exists, our goal is to show its existence. Hence our approach can be regarded as 
working in a kind of Lagrangian coordinates [16]; assuming that the solution exists we rewrite 
the system in the new variables ip u +u induced by u through (2.7). Then in the new coordinates 
we hope to be able to apply a fixed point method to show the existence of a solution. Since (u, w) 
are perturbations that are assumed to be small, we can assume that ||u||vk 2 is small enough that 
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the assumptions of Lemma 3 are satisfied. Hence the solution gives a well defined transformation 
x = ip u+uo (z), that we denote for simplicity by ip. If we show in addition its uniqueness in a class 
of small perturbations, then denoting = we have il>(Q u ) = ^ an d we come back to the 
original coordinates where our solution solves (2.5). Rewriting the system (2.5) in coordinates z 
yields 

u • V ' Z V + (V p o ip u ) d Zl u — fiA z u — (/i + z/)V 2 div 2 u + 7V 2 w = F(u, w), 
(V p o fa) d Zl w + div 2 u = G(u, w), 

n ■ 2[xD z (u) ■ r k + f(u ■ T k )\r = B k - n ■ 2/iR(u, D) • r k , 
n-u| r = 0, w\r m =w in . 

The functions F and G involve commutators from the change of variables. More precisely, 

F(u, w) = F(u,w)-u- R(V, V) - V p R{u, d Xl ) + fiR(u, A) + (/i + u)R(u, Vdiv) - -fR(w, V) 

(2.12) 

and 

G(u,w)=G(u,w)-R(u,dxv). (2.13) 

Here the first variable in the commutator R(-,-) denotes a function and the second is a differential 
operator. For example, R(w, V) := V x w — V 2 w and its i-th coordinate reads 

rW( w , v) = [9 2 ^(0« - 1) + X>^?]- 

We shall not give here precise formulas for the other commutators. Instead, we are now ready 
to give some heuristic arguments that will show what regularity we expect from the change of 
variables 0. To this end note that the commutators of the operators of order k depends on the 
derivatives of up to order k. More precisely, commutators of order one contain only compo- 
nents of the form 

while second-order commutators contain the terms 

V^-(V^) 2 , V z u-V 2 J 

and the terms of lower order. Hence in order to find the estimates on 1 1 F(u, w) \ \ Lp and 1 1 G(u, w) \ \ w i 
we need 

V x (f) G Loo, G L p , 

what will be satisfied provided that e W£ due to the imbedding W p l G (recall that p > 3). 
We should also note that, for simplicity of notation, the functions F(u, w) and G(u, w) in (2.12) 
and (2.13) denote exactly the same quantities as before. Hence we should keep in mind that they 
contain differential operators and so now they also contain some commutators that we will have 
to control to repeat the estimate (2.6). 

In order to show such estimate we should have a closer look at the derivatives of 0. With our 
construction of ip — 0~\ it is easier to consider first the derivatives of ip that are given by (2.10) 
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with u = ^yf 1 where u is the solution to (2.5), not to be confused with u = yp from the proof 
of Lemma 3 where u was an arbitrary function. In order to find the bounds on the commutators 
we will need the following smallness results for our change of variables: 

Lemma 4. We have 

^ H^T " ^ + 2. H^-ll^ 1 ^ ' (2 - 14) 

\\D 2 M Lp <E, (2.15) 

Enl-- 1 ik+Eiii-ik^. c-io 

\\D 2 J\\ Lp <E, (2.17) 

where E is sufficiently small number comparing to ji, depending on norms of perturbations 
measured by D defined in (1.11). 

Proof. The core of the proof is in the imbedding C L^. We start with (2.14). We 
estimate L p norm of the entries of E (2.10). These results quite directly from the form of E, but 
needs certain attention as E depends on ip implicitly. The entries of E without integrals will be 
small provided that ip is bounded what obviously holds true. For the entries involving integrals 
we change the order of integration and derivative obtaining 

d Zi 4> {j) = d Zi [ 1 u(^(s,z))ds= [ Z \w x u^(ij(s,z))-d z ^( S ,z)}d S . 
Jo Jo 

By Jensen inequality we have 

f 1 \dx k u^^(s,z))\ P \d z ^ (k \s,z)]\ P ds < Wr^VMl^ P \d Zi ^ k \s,z)]\ P ds. 
Jo Jo 

Integrating the last inequality over f2 we get 

ll^ylk < C(l + J2 II^IUJIIV^H^. (2.18) 

k,l 

The smallness of u in W£ and the imbedding W p C Loo gives (2.14). 

To show (2.15) we differentiate the entries of E, let us focus on entries with integrals. We 
have (we omit the sum over k): 

d zi d z ^ (j) = 9 ZI [ Zl [d Xk u^^(s,z))d Zi ^ k \s,z)}ds = 
Jo 

= [ 1 8 Zl [d Xk u^(tP(s, z))]d z ^ k \s, z)ds+ / 1 d Xk u^(s, z)d Zl d Zk ^ k \s, z) ds =: h + h- 
Jo Jo 
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Again by Jensen inequality, 



(2.19) 



and 

r>Zl 



w < nv^r Loo H \d zl d Zk ^ k \ s ,z)\ p ds. 

Jo 

Integrating the above inequalities over f2 we arrive at 

||V 2 ^||l p < C(fi)[||V^|| Loo ||V 2 ^|U p + ||V 2 «|U P ], (2.20) 

where the term HV^II! from (2.19) has been put into the constant. Like in the previous 
estimate, the imbedding C and the smallness of u in W 2 yield (2.15). 

To show (2.16) note that the smallness of E given by (2.14) combined with the imbedding 
Wp C Loo implies that detDtp > c > 0, hence Dip is invertible and we have 

D(f) = Dip- 1 = Id + E, 

where the elements of E can be explicitly computed in terms of E. The smallness of u together 
with the fact that is and algebra implies smallness in L p of the entries of E, which gives 
(2.16). Finally (2.17) is obtained by taking derivatives if D(j) similarly to (2.15). □ 

Now we are ready to show the basic estimate on the r.h.s. of (2.1 1): 

Lemma 5. Let F and G be defined by (2.12) and (2.13). Then we have 

\\F(u,w)\\ Lp + \\G(u,w)\\ w i < C[(\\u\\ W 2 + \\w\\ w i) 2 + D ] + E (\\u\\ W 2 + \\w\\ w i). (2.21) 

Proof. The bound on \\F(u, w) \\ Lp + \\G(u, w) \\ w i results from (2.6) in Lemma 2 (there are 
also some commutators since F and G involve differential operators, but these can be estimated 
as follows). We briefly justify the bounds on the commutators in F. To start with, the first order 
commutators contain the given function V, the functions u and w and the derivatives of 0, but 
only of the form 

4>%i±j and 0-1. (2.22) 

Hence applying Lemma 4 we get 

\\u- R(V, V)\\ Lp + \\V p R(u,d Xl )\\ Lp + \\R(w, V)\\ Lp <E(\\u\\ Wp2 + \\w\\ w i). (2.23) 

The second order commutators contain the second order derivatives of and the first order deriva- 
tives only of the form (2.22). Hence the application of Lemma (4) yields 

\\R(u, A)\\ Lp + \\R(u,Vdiv)\\ Lp < E \\u\\ W 2 (2.24) 

and we conclude the bound on In order to estimate \\G(u, w)\\ w i we differentiate the 

commutator 

i2(«,div) = 5>«(0« - 1) +5>g ) 0g ) (:r), 
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what yields 

o Zk R(u, div) = HIM - 1) + m ? E € 



Applying again Lemma 4 we get 

||i2(u,div)|| % i < £7||^||w^ (2.25) 

and the proof is complete. □ 

From now on we focus on the system (2.11) instead of (2.5). It is of crucial importance for us 
that we can solve (2. 1 1) in the domain H, what results from our choice of the transformation ip. 
Now we are in a position to define the operator 

T : W*{tt) x Wffi) ->• W*(tt) x Wffi), (2.26) 

to which we want to apply the Banach fixed point theorem in order to solve the system (2.1 1). 
Namely, we set (u, w) = T(u, w) if 

u ■ V Z V + (V p o ipu+u ) d Zl u - fiA z u - (/i + v)V z div z u + -fV z w = F(u,w), 
((V p o ipu +ua )) d Zl w + div z u = G(u, w), 

(2.27) 

n ■ 2/iD z (w) • r k + f(u ■ T k )\ r = B k , 
n-u|r = 0, w\r m =w in . 

The point is that the term d Xl w + u ■ Vw is replaced by ((V p o + u^)d zi w, and for this term 
we find a bound in W^, what is necessary to show the contraction property of T. From now on 
for simplicity we will write tp a instead of ipu +Uo . 



3 A priori bounds and solution of the linear system 

In this section we deal with the linear system: 

u ■ V ' Z V + (V p o ip u ) d zi u — fiA z u — (/i + v) V z div z u + w — F, 
((V p oi, u ))d Zl w + div z u = G, 

n ■ 2fjD z (u) ■ r k + f(u ■ T fc )| r = B k , 
n ■ u\ r = 0, w\ Tin = w in , 

with given functions F, G, B k , u. Note that we take the superposition of V p and ^« to obtain 
V p o%jj u for the original system and hence V p when we pass to the original system of coordinates. 
The same remark concerns the function F, but it does not change anything in the computations. 
We need to solve the system (3.1) to show that T is well defined by (2.27). To this end we need 
the appropriate estimates that we show in the first part of this section. In the second part the 
linear system is solved. 
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3.1 A priori bounds 

In this section we show the estimate in W 2 x for the solution of the linear system (3.1) in 
the maximal regularity regime. The first step is the energy estimate. It is given by the following 

Lemma 6. Let (u,w) be a solution to the system (3.1) with given (F,G, B,u) G V* x L 2 x 
L 2 (T) x Wp, where u is small enough to assure 

V p + u w > 9 3 > 0, (3.2) 

for some 9 3 > 0. Assume that the friction f is large enough on Y in and the viscosity /x and the 
friction on r satisfy (1.12). Then 

Nlw* + IMUoo(l 2 ) < C [\\F\\v + ||G||l 2 + ||5|| L2( r) + lkm||L 2 (r in )], (3.3) 

where 

V = {v e W£(Q.) : v ■ n\ T = 0} (3.4) 

and V* is the dual space ofV. 

Proof. We start with two basic observations. First, since V£ = 0, by (2.7) we have 

d Zl (V P o ifo) = ±;[u^d X2 V p + u^d X3 V p }. 

Now recall that the constant 9 in (1.5) is independent of the smallness of perturbation. Hence we 
can assume the ||u|| W 2 is small compared to 9 and by the imbedding C Loo we have 

\\d Zl (V p o^)\\ K <E(u,V p ). (3.5) 

We keep in mind that (1.12) and (1.6) hold, what implies that (3.5) will be controlled by the 
viscosity, more precisely even by a constant which decreases with increasing viscosity due to 
(1.6). In the remaining of this section we will write V p instead of V p o The fact that we 
consider the superposition does not influence the computations as we have (3.5). We apply the 
identities 

/ V p d Zl \u\ 2 dx = I I V p \u\ 2 n (1) da - [ \u\ 2 d Zl V p dx (3.6) 



(3.7) 



and 

/ (— /jlAu — (u + /i)V divu) • v dx — / 2/iD(-u) : Vf + u divu div v dx— 
Jn Jn 

/ n ■ [2fiD(u)] ■ vda — / n ■ [u(divu)ld] ■ v da. 

Now we multiply (3.1)i by u and integrate. Using the above identities, with application of the 
boundary conditions (3.1) 3)4 we arrive at 

f {2fjD(u) : D(u) + v\div u\ 2 } dx + j (f - X -^-)\u\ 2 da + j + V p )\u\ 2 da 

—7 / w div udx+ / (u • VV) ■ u dx — (3.8) 
Jn 

= / F-udx+ / {Bi(u ■ n) + B 2 (u ■ r 2 )} da + / \u\ 2 d Zl V p dx. 
Jn Jr Jn 
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Note that by (3.5) we have 

| / \u\ 2 d Zl V p dx\ < E(u,VV p )\\u\\i 2 . 
Jn 

The other terms on the r.h.s. are all 'good' terms. The T out term on the l.h.s. is nonnegative and 
the T in term will be positive for / large enough. To deal with the term f Q w div u dx we apply 
the continuity equation to express div w obtaining: 

/ wdiv u dx = — I Gw dx — - j w 2 d Zl (V p + u^) dx 
Jn Jn 2 J n 

+ 1 [ W \V P + u^)da-\ I w 2 n {V p + u^)da. 
& /p z /p. 

° 1 out ° J- tn 

By (3.2) the integral over T out will be nonnegative, hence we have 

/ wdW U dx< ||G|| L J| W || L2 + ||9 2l (\/ p + ^)|| Loo || W ||i 2 +C|| Wm ||i 2(rm) . (3.9) 
Jn 

To derive the W\ - norm of u we apply the Korn inequality [23, 27]: 

/ [2/xD(u) : D(m) + u\d\Yu\ 2 }dx+ / / (u • r) 2 rfa > C K \\u\\ 2 wl , (3.10) 

where = C k (/jl, u, f, Q) and Cx is increasing with /i. A sketch of the proof of (3.10) one can 
find in the Appendix, note that in (3.10) we use only information at T in , a part of the boundary, 
but still it is sufficient to control the whole norm of W\. 
Combining (3.8), (3.9) and (3.10) we get 

C K \\u\\ 2 w i < \\G\\ L2 \\w\\ L2 + C\\w in \\ 2 L2{r . n) 

+ [||F||^ + || J B|| L2(r) ]|| M || H , 2l + J E(M,V\/ p )[||^||| 2 + || M ||| 2 ]- / (u-VV)-udx. (3 - 1!) 

Jn 

We have to deal with the last term on the r.h.s. It is impossible to show it has a good sign, hence 
the only way is to estimate it directly with 

| / (u-VV)-udx\ < C 2 P \\VV\\ L J\u\\ 2 j, (3.12) 

J Si 



where Cp is the constant from the Poincare inequality in V. Inserting the above to (3.1 1) we get 

(3.13) 



{C K - C\ (IIVVIU. + E(u, VV P )) \\u\\ 2 < 



< E |MI! 2 + C(\\G\\ L2 + ||«; in || L2( r in) ) + (||F|| V . + \\B\\ L2{T) )\\u\\ w i, 

where we recall that E is a small constant and C is a data-dependent constant, not necessarily 
small. Now, Ck is increasing with fi, while Cp does not depend on ji. Moreover, (1.6) implies 
that E(u, W p ) will be decreasing when ji increases. Finally, (1.8) implies that for small values 
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of the friction / on T it is enough to assume that the viscosity is large only compared to /|r 
to control E(u, W p ). We conclude that the constant on the l.h.s. of (3.13) will be positive 
provided that /j, satisfies (1.12). 

Here it is a good point to emphasize the necessity of sufficient magnitude of the viscosity 
coefficient for large / at r . This assumption is somehow natural, although in the case V p = 
const it is not required [19]. We have to control (3.12), and largeness of dissipation may only 
give us this chance. Note that for the Dirichlet boundary condition [9], although the constant 
flow is considered, such assumption is required, too. 

To complete the proof of (3.3) we find a bound on HHIw^)- To this end we refer to the 
next subsection, where we solve the linear system. Notice that w = S(G — divu), where S is 
defined in (3.35), and so 

HlLco^^CdlGlk + IN^). (3-14) 
Combining (3.13) and (3.14) we conclude (3.3). □ 

In the next step we show higher bound on the vorticity of the velocity. To this end we take 
the vorticity of (3.1)i. Denoting a = rot u we get 

-fiAa = rot [F - V p d Zl u - u ■ VV] in SI, 

a-r 2 = (2xi - f -)u ■ ti + ^ on T, 

a .T 1 = (£-2 X 2)u-T 2 -f on I\ ^ 

div a = on T. 

The boundary conditions (3.15)2,3 are derived from differentiation of (3.1)4 in tangential di- 
rections and application of (3.1) 3 , see [15], [19]. The above system gives the estimate ([27], 
Theorem 10.3 with \x = 0): 

IMIw? <C[\\F\\ Lp + H^Hw^lMlw^ + \\u\\ w i-i/ P{r) + \\B\\ w i-i /P(r) } 

< c [\\ F h P + \\B\\ w i-i /p+p{r) + ||u||wi]. 
Applying the interpolation inequality (5.7) to ||m||w/i and the energy estimate we get 

IMIwi < C( e ) [\\ f \\l p + \\G\\ w i + \\ B \\ w i-i/ P{r) + \\win\\wi(r in )] + 4 u \\w^ (3-16) 

for any e > 0. Now consider the Helmholtz decomposition of the velocity 

u = V(f> + A, (3.17) 

where |^| r = and div A = 0. We see that the field A satisfies the system 

rot A = a in f2, 

divA = in SI, (3.18) 
A ■ n = on T. 

This is the classical rot-div system and from [23] we have ||^4||vk p 2 < C \\a\\ w i, what by (3.16) 
can be rewritten as 

\\A\\ W 2 < C(e) [\\F\\ Lp + \\G\\ w i + \\B\\ w i-i /P(r) + ||^ n || Wp i (ri „)] + e\\u\\ W 2 (3.19) 
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for any e > 0. Now we substitute the Helmholtz decomposition to (3.1)i. We get 

V[-(z/ + 2 / u)A0 + 7w] = 

(3.20) 

F - V p d zl A + fxAA + (v + /x) V div A - A ■ W - V p d Zl V<p - V0 • W =: F. 
Since A0 = div u, we can write 

-(z/ + 2/i)divw + 7w = H. (3.21) 

We underline that we are now at the level of a priori estimates and (3.21) should be treated as the 
definition of H . In fact we can think of H as a kind of effective viscous flux like in the theory of 
weak solutions to compressible Navier-Stokes equations ([7], [12]). Now (3.20) can be rewritten 
as VH = F. Combining the last equation with (3.1) 2 we arrive at 

^w + V p d Zl w = H, (3.22) 

where 7 = — V and 

H=^— + G. (3.23) 

The equation (3.22) makes it possible to estimate IIHIwp- an d H^ziHIvKp 1 m terms of H. Next 
we can find the bound on H using interpolation and the energy estimate. The first step is in the 
following lemma: 

Lemma 7. Let w solve (3.22) with H e andw\r in = w in . Then 

\\ w \\w^ + \\d Zl w\\w^ < C [\\H\\ w i + \\w in \\ w i {rin) }. (3.24) 

Proof. To estimate ||iu||l p we multiply (3.22) by \w\ p ~ 2 w and integrate. Using the boundary 
conditions we get 

tIIHI? + - f \w\ p v p d(r< 



< HfflklkllC + - / \w\ p V p da+ (\\d zl V p \\ Loo + \\u\\ wL ) Ml' 1 . 



The boundary term on the l.h.s. is positive and the constant in the last term on the r.h.s. is small 
(note that we take only the z 1 derivative of V p ). Hence the above implies 

lhl|L p <C[||//|| ip + ||^ n || M r in) ]. (3.25) 

In order to find a bound on d Zi w we differentiate (3.22) with respect to z%. If we assume that 
w G Wp then (3.22) implies d Zl w G since Wp is an algebra. Thus we differentiate (3.22) 
with respect to Zi, multiply by \d Zi w\ p ~ 2 d Zi w and integrate. We have 
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f (y p + #) \d Zi w\ P ~ 2 d Zi wd Zl w = - [ {V p + u {1) )d Zl \d Zt w\ P dx = 
Jn P Jn 

{V p + u^)\d Zi w\ P n {1) da- I \d Zi w\ p d Zl (V p + u {l 
Jt Jn 



(3.26) 



The last term on the r.h.s. can be estimated by E \\d Zi w\\ Lp since we take only z\ derivative of 
V p . The boundary term vanishes on T and Y out part will be nonnegative. Hence (3.26) implies 

l\\d Zi w\\l p < \\H\\w} \\d Zi w\\ P C + ~ f \d Zi w in \ p (V p + u^) da + E\\d Zi w\\l p - 0-27) 

For % = 2, 3 we use the fact that w in e W p l (Y in ) and conclude that 

R.Hk < C [\\H\\ W , + \\w in \\ n{Tin) ]. (3.28) 

To apply this method to d Zl w we need some knowledge on d zi w in \r in . To this end we can use 
(3.22), which, since V p + > 0, can be rewritten on Y in as 

H - ^w in 

C Zl W in ~ yp + • 

Hence \\d Zl w\\ Lp(Tin) < C[\\H\\ Lp(Tin) + \\w in \\ w i {Tin) ], and (3.27) implies (3.28) also fori = 1. 
From (3.25) and (3.28) we conclude 

\\w\\ w i < C [\\H\\ Lp + ||^ n || Lp (r in )]. 

The bound on ||9 Xi k;||v^ results simply from the identity (3.22) and the fact that W p x is an algebra. 
The proof of (3.24) is complete. □ 

Now we need to find the bound on \\H\\ w i, but this is straightforward. Interpolation inequal- 
ity (5.7) yields 

\\h\\ Lp <s\\vh\\ Lp + c(S)\\h\\ L2 , 

for any 5 > 0. To estimate \\H \\ L2 we use the fact that H — + 2fi)div u + and the energy 
estimate (3.3). To find the bound on ||Vif||x, p we use (3.19), (3.20), then (5.7) to estimate the 
term IHInn and finally (3.3). We obtain 

\\H\\ W } < C(S) [\\F\\ Lp + \\G\\ w i + \\B\\ w i-i /P(r) + ll^inllwri^)] + 5\\u\\ W 2. (3.29) 

We are now one short step from the main result of this section. It is given by the following 

Theorem 2. Let (u, w) be a solution to (3.1) with (F, G, B k ,w in , u) G L p x W p l x Wp~ 1/p (T) x 
Wp(T in ) x Wp such that \\u\\ W 2 is small enough and the friction f is large enough on T in . 
Assume also that the viscosity and the friction on T satisfy (1.12). Then 

IMIw? + \\ w \\w} + \\d gl w\\w} < C [\\F\\ Lp + \\G\\ w i + \\B k \\ w i-i /P{rin) + \\w\\ w i {rin) }. (3.30) 
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Proof. To close the estimate (3.30) it remains to find the bound on ||m||w|- To this end notice 
that in particular u satisfies the Lame system: 

—/j,Au — (u + fi) Vdiv u = F — 7Vu> — V p d Zl u — u ■ VV in Q, 

n ■ 2/iD(w) • ^ + f u ■ n = B h i = l,2 on T, (3.31) 

n ■ u = on T. 

Lemma 1 1 applied to the above system yields 

IMIw| < C[\\F\\ Lp + \\w\\w$ + \\ B \\ w $-v*(r) + hWw*]- 

Applying the interpolation inequality to the term and then the energy estimate (3.3) we 

get 

\\u\\wl < C [\\F\\ Lp + \\G\\ w i + \\w\\ w i + \\B\\ w i-i/ P{r) + \\w in \\ L2(Tin) \. (3.32) 

Combining this estimate with (3.24) and (3.29) with appropriate 8 we conclude (3.30). □ 

Note that Theorem 2 gives more than we need to solve (3.1), namely the bound in Wp of 
d Zl w. We will use this result to show the contraction property for the operator T in Section 4. 

3.2 Solution of the linear system 

With the estimates that we obtained we are ready to solve the system (3.1). First we define the 
weak solution and show its existence. Next, applying the estimate (3.30) we show its regularity 
under the appropriate regularity of the data. 

3.2.1 Weak solution 

By the weak solution to (3.1) we mean a couple (u, w) E V x L (X> (L 2 ) such that 

{v ■ V p d Zl u + u ■ VV) + 2/iD(u) : Vi> + udivudivv — 'ywdivv} dx 

+ / f( u ■ T i) (v ■ Ti) da = / F ■ vdx + / B^v ■ t { ) da (3.33) 
Jr Jn Jr 

is satisfied V v G V and (3.1) 2 is satisfied in V'{Q), i.e. for all G C°°{Q), <f)\r out = 0: 

- / (V p +u^)wd Zl <f)dx- [ d Zl (V p +u^)w(Pdx = [ (j)(G-divu)dx+ [ {V p +u^)w in (t)da. 
Jn Jn Jn Jr in 

(3.34) 

To find the weak solution we apply the Galerkin method. Hence we introduce an orthonormal 
basis of u k C V and finite dimensional spaces V N = {J2iLi a i u i '■ <^i e M.} C V . We look 

N 

for the approximations of the velocity of the form u N = cfuj. Taking into account the 
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continuity equation we have to define the approximations of the density in an appropriate way. 
Namely, we set w N = S(G N — div u N ), where S : L 2 (Vi) — > L OQ (L 2 ) is defined as 

w = S (v) ~ ( a" + *»>*.— - D-(n). (3 . 35) 
{ w = w in on T in 

and satisfies the estimate 

l|5(t;)|| L(w(La )<C7[||Ti; in || LaC r iB ) + |H| La ]. (3.36) 
The construction of S is quite straightforward. For a continuous v we set 

S(v)(z) = w in (0,z 2 ,z 3 ) + J vp + _ w (s,z 2 ,z 3 )ds. (3.37) 

Next we show directly the estimate (3.36), which enables us to extend S on L 2 (VL) using a 
standard density argument. 

Now we proceed with the Galerkin scheme. Taking F = F N , u = u N = Yli c i v — 
Uk, k — 1 . . . N and w = w N = S(G N — div u N ) in (3.33), where F N and G N are orthogonal 
projections of F and G on V N , we arrive at a system of iV equations 

B N (u N ,u k ) = 0, k = l...N, (3.38) 

where B N : V N ->■ \/ w is defined as 

B N (£ N , v N ) = fn{v N V p d Zl £ N + £ N ■ VV + 2/iD(e iV ) : Vu^ + div ^ div u^} 
-7 / n - div^) div / rfx + / r [/ • r,) - 5,] (v N ■ Tj ) da — J n F N ■ v N dx. 

(3.39) 

Now, if u N satisfies (3.38) for A; = 1...N and w N = S(G N - divw^), then a pair (u N , w N ) 
satisfies (3.33) - (3.34) for (v,<j>) G (V N x (j)\r out = 0. We will call such a pair 

an approximate solution to (3.33) - (3.34). To solve the system (3.38) we apply the following 
well-known result (the proof can be found in [25]): 

Lemma 8. Let X be a finite dimensional Hilbert space and let P : X — >■ X be a continuous 
operator satisfying 

3M>0: (P(O,O>0 for ||£|| = M. (3.40) 
Then 3f* : ||f || < M anJ P(^) = 0. 
We define P N : V N ->■ \/ w as 

^(^) = E 5 ^^H for {"6^. (3.41) 

In order to apply Lemma 8 we show that (P(£ N ), £ N ) > on some sphere in V N with radius 
dependent on the norms of the data. To this end we follow the proof of the energy estimate for 
(3.1). This is in fact standard approach in the Galerkin method: the energy estimate combined 
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with Lemma 8 gives existence of the approximate solutions, hence we skip the details here. 
Except from the existence of the approximate solution u N , Lemma 8 gives the estimate 



u N \\ w i < C(DATA), 



which combined with (3.36) gives 



u \\ w i + || w 



N 



Lao{ L 2) < C(DATA). 



Thus 



u 



N 



u in W 2 and w 



N 



* w in ^(Ls) 



for some (u, w) E W\ x L 00 (L 2 ). We easily to verify that (u, w) is a weak solution. First, passing 
to the limit in (3.33) for (u N , w N ) we see that u satisfies (3.33) with w. On the other hand, taking 
the limit in (3.34) we verify that w = S(G — div u). We conclude that (u, w) satisfies (3.33) - 
(3.34), thus we have the weak solution. To show the boundary condition on the density we can 
rewrite the r.h.s of (3.35) as 



and, treating x\ as a 'time' variable, adapt Di Perna - Lions theory of transport equation ([6]) that 
implies the uniqueness of solution to (3.42) in the class L OQ (L 2 ) (note that this is the reason we 
work with weak solutions with the density in L 00 (L 2 ) instead of usual L 2 ). This completes the 
proof of existence of weak solution. 

3.2.2 Strong solution 

The following result gives strong solution to the linear system (3.1) for the data of appropriate 
regularity. 

Theorem 3. Let (F,G,B k ,u) e (L p x x Wp~ 1/p (T) x W%) with \\ u\\ W 2 small enough. 
Assume further that f is large enough on Y in and jj, fulfills (1.12). Then there exist a unique 
solution (u, w) G W£ x W p to the system (3.1) and the estimate (3.30) holds. 

Proof. To show appropriate regularity of the weak solution for the regular data it is enough 
to apply the estimate (3.30) provided that we handle the singularities of the boundary at the 
junctions of the wall r with inlet T in and outlet T out . To this end we apply the result on the 
elliptic regularity of the Lame system with slip boundary conditions, Lemma 1 1 in the Appendix. 
Notice that we can apply this method since we work in the fixed domain VI due to the appropriate 
choice of the change of variables. Otherwise we would end up in a free boundary problem and 
the solution of the linear system would be highly nontrivial. □ 




(3.42) 
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4 Contraction 



In this section we show the contraction property for the operator T defined in (2.27). However, 
first of all we notice that Theorem 3 gives the solution of the linear system (3.1) provided that 
||«|| W 2 is small enough, and so this constraint must hold if we want to have T(u, w) well defined. 
We start this section with showing a stronger result, namely that T : B R — > B R for some ball 
B R e Wp x Wp with R depending on the data. This property will also be needed to show the 
contraction. Eventually, we prove Theorem 1, the main result of the paper. 

Lemma 9. There exists R > depending on the size of the data measured by Do (1-11) such that, 
provided the data is small enough, T(B R ) C B R , where B R is a ball of radius R in Wp x W p l . 

Proof The estimates (2.21) and (3.30) imply that 

\\T(u,w)\\ W 2 xW i <C\\(u,w)\\ 2 W 2 xW i+5, 
where 5 < provided that the data is small enough. For such 5 we have 

II (u,w) \\w*xwi < 25 \\T(u,w)\\ W 2 xW i < 25. □ 

Now we show the contraction property for T. To this end consider = T(ui,wi), 

(u 2 ,w 2 ) = T(u 2 ,u> 2 ). By the definition of T, the difference (u,w) := (ui — u 2 ,w 1 — w 2 ) 
satisfies the system 

u ■ V Z V + V p d Zl u - fiA z u — + v)\7 z div z u + 7V \ w = F(ui,wi) - F(u 2) w 2 ), 
(V p + u { 2 1] ) d zi w + div z u = G(ui,«)i) -G(u 2 ,w 2 ) + (ui -u 2 )d zi w u 

n ■ 2/iD z (u) ■ r k + f(u ■ T k )\ r = n ■ 2/i[R(u 2 ,'D) - R(ui, D)] • r k , 
n-u\ r = 0, w\r m = 0- 

Hence in order to show the contraction principle for T we can apply (3.30) provided that we have 
good bounds on the r.h.s. of (4.1). A result we need is given by the following 

Lemma 10. We have 

\\F(ui, wi) - F(u 2 ,w 2 )\\ Lp + \\G(u!,ivi) - G(u 2 ,w 2 )\\ w i 
+ ||(«i - u^d^WxWwi + \\n ■ 2(x[R(u 2 ,~D) - R(u 1: D)] • Tfcll^i-i/p^ (4.2) 
< E [Wm - u 2 \\ W 2 + \\wi - «; 2 ||w-i]. 

Proof Denote (u i: Wi) by (u i: Wi) for simplicity. We have 

F(ui,wi) - F(u 2 ,w 2 ) = Mi • V x u x - u 2 ■ V x u 2 + (ui - u 2 ) ■ V x u + u ■ V x (ui - u 2 ) 
+5ti'(wi)'V x Wi - 5ti'{w 2 )V x w 2 + (w l - w 2 )F + R, 

where 

5n'(wi) := -n\wi + 1) - n'{wi) 
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and R denotes all the differences between corresponding commutators in F. We can afford using 
such abbreviation and estimate R without any additional computation if we just notice that the 
commutators are linear with respect to the functions and hence if certain estimate in terms of the 
function holds for a commutator, then the same estimate in terms of the difference hold for the 
difference of the commutators, for example 

\R(u,d Xl )\ < E\\u\\ w i \R(ui,d Xl ) - R(u 2 ,d Xl )\ < E \\m -u 2 \\ w i. 
Applying this reasoning to all the commutators in F we conclude that 

\\R\\l p < E(4>)(\\ Ul -u 2 \\ W 2 + \\ Wl - w 2 \\ w i). (4.3) 
We estimate the remaining parts. Obviously we have 

||(ui - u 2 ) ■ V x u \\ Lp + ||m • V x (mi - u 2 )\\ Lp < £(w o ,0)IK - u 2 \\ W 2, (4.4) 

where E(-, •) depends also on since we have commutators as the gradients are w.r.t. x. A little 
bit closer examination shows that the i-th coordinate 

(ui ■ V x ui - u 2 ■ V x u 2 ) W = (ui ■ V x (ui - u 2 ) + (ui - u 2 ) ■ V x u 2 ) W = 

= [fa - «2)°' ) ^t4 i) + r(4\ d Xj )] + J2 («i - + ^((«i - u 2) (i \ d Xj )} , 

3 3 

and so by a direct computation we get 

IK • - u 2 ■ V x u 2 \\ Lp < E{\\ui\\ W 2,<f>)\\ui - u 2 \\ W 2. (4.5) 
It remains to estimate 

Sir' '(wi)V >i - 5it' (w 2 )V x w 2 = 5tt'(wi)V x (wi - w 2 ) + [5n'(wi) - 5n'(w 2 )] V x w 2 = 

= 8tt'(wi)[V z (w 1 -W2) + R(w 1 -W2,V)] + |fc'(t«i) - Sir' (w 2 )][V z w 2 + R(w 2 , V)]. 
It follows easily that 

||57r'(wi)V x wi - 6n'(w 2 )V x w 2 \\ Lp < EQlimWwi, \\w 2 \\ w i, <j>)\\wi -w 2 \\ w i. (4.6) 
Combining (4.3), (4.4) , (4.5) and (4.6) we conclude 

||F(wi,wi) - F(u 2 ,w 2 )\\ Lp < E(\\wi\\ w i, ||^||w2,0)[||wi -u 2 \\ W 2 + \\w! -w 2 \\ w i]. (4.7) 
Now we estimate the difference in G. We have 

G(u!, wi) - G(u 2 ,w 2 ) = 
= (w 2 — wi)divu + wi div(-ui — u 2 ) + d\YU 2 {w\ — w 2 ) + R(u 2 , div) — R(u±, div). 
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The first term 

\\(w 2 - Wi)divwo||w p i < ||V(w 2 - ^i)divM ||L p + ||(w 2 - Wi)V 2 m ||l p < 

< ||divuo||Lj|wi -w 2 \\ w i + ll^o || w-2 ||^i -w 2 \\ Loo < C|ko||w2|l w i - w 2 \\w^- 
The second 

\\w 1 div( y u 1 -u 2 )\\w} < ||(Vwi)div(Mi-M 2 )||L p + lki V 2 (mi-m 2 )||l p < C \\w 1 \\ w i\\u 1 -u 2 \\w^ 
Similarly we show 

||divM 2 (wi - w 2 )\\ w i < C|M|w| IK -w 2 \\ w i. 

Now we have to estimate the difference of the commutators in W*. It turns out to be straightfor- 
ward as we have 

R( Ul , div) - R(u 2 , div) = £ d Zi (ui - u 2 f>0 + £ d Zj (uU - uY ] 

i i^j 

and so, identically as in (2.25), we show that 

||i2(«i,div) - R(u 2: div)\\ w i < E{4>)\\{ Ul -u 2 )\\ w ,. 
Combining the above results we conclude 

||G(«i, w ± ) -G(u 2 , w 2 ) || w i < [E((f)) + C (\\w4wi + ||u 2 ||w*)] (IK - u 2\\w* + \\wi-w 2 \\ w i). 

(4.8) 

To treat the last term of the r.h.s. of (4.1) 2 we observe that 

|| («i - U2)d Zl wi\\ w i < CWdztWxWwiWv,! -u 2 \\ W 2. (4.9) 

It remains to estimate the boundary terms n ■ [R(u 1: D) — R(u 2 , D)] • r k . To this end it is enough 
to notice that 

{R(u l7 D) - R(u 2 , D)}ij = R(( Ul - u 2 f\ d Xj ) - R(( Ul - u 2 )^\d Xi ) =: R hr 
Applying the trace theorem and repeating the proof of (4.8) we can show that 

ll^ijllv^i-i/^r) < C Il-Rijllwp 1 < E{(f>) \\ui -u 2 \\ W 2, 

and so 

\\n ■ [i?K,D) - R(u 2 ,T>)] ■ T k \\ w i-i /P(r) < E(4>) |K -u 2 \\ W 2. (4.10) 

Combining (4.7), (4.8), (4.9) and (4.10) we conclude (4.2). □ 

Proof of Theorem 1. With the results of the previous section we can apply the Banach 
fixed point theorem to the operator T. It gives existence of a unique fixed point within the 
ball B(0, R) e Wl x W*. By Lemma 9 we have R = R(D ) where D is defined in (1.11). By 
the definition of T, the fixed point (u, w) solves the system (2. 11). 

Now we recall Section 2 and conclude that the original coordinate system is x = ip u (z), and 
in the x variable our solution satisfies the system (2.5). It follows that v = v + u + u$ and 
p = w + 1 solves (1.1) and the estimate (1.13) holds. Theorem 1 is proved. 
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5 Appendix 



Proof of Lemma 1. As explained in Section 1, we assume the pressure of the form TL = u(f) x±. 
Then on each x\ - cut of f2 (i.e. on each set f2 x {xi}) V p can be found as a solution to the 
elliptic problem (1.7) which we recall here: 



/lAv = U Xl = ui(f) < in f2 
^ + /« = on an . 



(5.1) 



Testing (1.7) with v^—v x v <o we get 



- / pt\Vv_\ 2 - / fvi = / u(f)v_ > 0. (5.2) 
Jn Jn Jn 

The last inequality results from nonpositivity of i>_ and the direction of the flow which implies 
u(f) < (recall the remark after (1.7)) Clearly (5.2) implies v > 0. We want to show sharp 
inequality. To this end consider v satisfying 

Av = u(f)/n<0 inn, v\ r = at dtt. 



We have v > 0, and, by the maximum principle applied to v — v we get inf^ v = inf r v. Assume 

dv 
dn 



that mi r v = v(x ) = for some x e T. Then, since / > and |^ < 0, we must have 



But since u(xo) = infn(u — u), by the Hopf Lemma we must have l9( -g~ t '- > (x ) < 0. The applica- 
tion of of Hopf lemma is possible since Q is a C 2 subset of IR 2 , hence v is a classical solution 
to (5.1), i.e. v e C 2 {Q) n C^O). We conclude that v > 6 > on T, hence v > 9 in Q. In 
particular, the linearity of (5.1) gives (1.6) where cu is continuous since cu is continuous (recall 
the discussion after the formulation of Lemma 1). This completes the proof. □ 

Lemma 11. (Lame system with slip boundary conditions). Let ji > 0, v + 2/i > 0, F e L P (Q) 
and B e W^p~ 1//p (r). T/zen ?/zere exz5?5 m G W^(Q) solving 

—fiAu—(u + fi)'Vdivu = F in Q, 

n-2fiD(u) ■ r + f u- r = B, 2 = 1,2 on V, (5.3) 
n • u — on Y. 

Moreover, the following estimate holds: 

\\u\\w><C[\\F\\ Lp + \\B\\ wr/P{T) l (5.4) 

Proo/ Under the assumptions on /i and z/ (5.3) is elliptic so we easily get a weak solution. The 
only problem we encounter showing regularity of the weak solution under appropriate regularity 
of the data are the singularities of the boundary on the junctions of r with Y in and T out . These 



25 



can be dealt with using symmetry arguments. But first we have to reformulate slightly the system 
(5.3). Having the weak solution, and the bound in W\ for the velocity, we are able to consider 
only the case / = and B = 0. The general case will be a consequence of this particular one. 
Assume that T in C {x 1 = 0}, then we define the operator E v as extending a vector field defined 
for {x : x\ > 0} on the whole space as 

where x = (— x-y, x 2 , x 3 ) and u(x) = [— u^(x), u 2 (x), u( 3 \x)]. Then we have 

AE v as (v) + VdwE v as (v) = Av + Vdivv, (5.6) 

and on the plane x\ = the extension E v as preserves the slip boundary conditions for / = and 
B = 0, since 

n ■ D(E: s (v)) ■ Ti\ Xl =o = d Ti E: s (v)?d Xl E: s (v)V = 

and n ■ E v as (v)\ Xl=0 = 0, d xl E v as (v)^ = for % = 2,3. Now we are allowed to localize the 
equations in a vicinity of T in obtaining a system in a smooth domain. Then the theory [2, 3] 
gives the full regularity of E^ s (v) in neighborhood of the junctions. Application of analogous 
antisymmetric extension on T out completes the proof. □ 

Lemma 12. (interpolation inequality): 
ForfeW^n), p>3: 

ll/IU P <6||V/|| Lp + C(e,p,fi)||/|| L2 . (5.7) 

Proof. The interpolation inequality in L p ([1], Theorem 2.11) and the imbedding W p C Loo 
for p > 3 yields 

ll/lk < c(p) Il/IIL 11/11^ < c(p)(ll/lk + llv/IUj e \\f\\]- d , 

what entails (5.7) after application of the Cauchy inequality. □ 

Lemma 13. Let Q be bounded with sufficiently smooth boundary and T part be an open regular 
subset of dQ. Then 

\\u\\ w i <C(||D(«)|| La + ||«| ) (5.8) 

foru e W£(n). 

Proof. The known result based on properties of the kernel of D(-) [1 1, 28] yields 

IMI^i < CiUD^H^ + CalHI^. (5.9) 
In order to obtain (5.8) we shall prove that 

\\u\\ L2 < ^tHV«||l 2 + M(||DH|| L2 + \\u\ rpart \\ L2(rpart) ). (5.10) 
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Compactness argument and features of D implies (5.10). We use that fact that the only solution 
to the system D(w*) =0, u\ Tpart = is u* = 0. To show it we use the fact that if D(-u) = then 




(5.11) 



Hence u is an affine map and it is easy to verify verify that if at least one of the coefficients 
A, B, C is nonzero then the rank of matrix of u is 2, hence Ker u is a line. On the other hand, 
T part is a two dimensional submanifold as it is an open, regular subset of dQ. But r part C Kera, 
hence we conclude that A = B = C = 0, hence di — 0. □ 
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